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ABSTRACT 



A promising mechanism to form intermediate-mass black holes (IMBHs) is the runaway merger 
in dense star clusters, where main-sequence stars collide and form a very massive star (VMS), 
which then collapses to a black hole. In this paper we study the effects of primordial mass 
segregation and the importance of the stellar initial mass function (IMF) on the runaway growth 
of VMSs using a dynamical Monte Carlo code for A^-body systems with N as high as 10® stars. 
Our code now includes an explicit treatment of all stellar collisions. We place special emphasis 
on the possibility of top-heavy IMFs, as observed in some very young massive clusters. We find 
that both primordial mass segregation and the shape of the IMF affect the rate of core collapse 
of star clusters and thus the time of the runaway. When we include primordial mass segregation 
we generally see a decrease in core collapse time (ice)- Although for smaller degrees of primordial 
mass segregation this decrease in tec is mostly due to the change in the density profile of the 
cluster, for highly mass-segregated (primordial) clusters, it is the increase in the average mass in 
the core which reduces central relaxation time decreasing tec- The final mass of the VMS formed 
is always close to ~ 10"'^ of the total cluster mass, in agreement with the previous studies and is 
reminiscent of the observed correlation between the central black hole mass and the bulge mass 
of the galaxies. As the degree of primordial mass segregation is increased, the mass of the VMS 
increases at most by a factor of 3. Flatter IMFs generally increase the average mass in the whole 
cluster, which increases t^c- For the range of IMFs investigated in this paper, this increase in t^c 
is to some degree balanced by stellar collisions, which accelerate core collapse. Thus there is no 
significant change in t^c for the somewhat flatter global IMFs observed in very young massive 
clusters. 
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Introduction 



1.1. IMBHs 



It is generally accepted that there exist two 
separate classes of black holes (BHs) defined 
by their mass ranges: stella r -mass BHs of mass 
^ 3 - - 20 MfT, dBoltonl [l972l: IWebster fc Murdin 
I1972I : ICasaresI l2007l ) that form through the col- 
lapse of massive stars, and supermassive BHs 
(SMBHs) of mass - 10® - 101°Mq that are 



found in the cent ers of m ost galaxies including 

the Milky Way (iKormen dv & Gebhardt 200l[ 

Merritt fc Ferraresel200l[lGhez et al.l2003HMiller fc Colbert 



2004I ). However, this leaves a gap in the mass 
range from ^ 10^ — 10^ M© which in recent 
years seems to have been filled by tentative ev- 
i dence for intermediat e-mass black holes (IMBHs ) 
(jPortegies Zwart et al. 1999: Grindlav et al.ll200ll ). 
This e vidence consists of dynamical measure- 
ments (ICebhardt et al.ll2002l: ICerssen et al.ll2002 : 
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van der Marel et al.l 120021 ). detection of ultralu- 
minous X-r ay sources (ULX) fou nd off-center 
in galaxies (jMiller fc Colberd 120041 ) , the anoma- 
lously high mass-to- light ratio ob seryed in the cen- 
ters of some globular clusters ([Miller fc Colbert 
I2OO4I) and the mass segregation quenching in the 
cores of globular clust ers in presence of an IMBH 
(|Pasquato et al.ll2009l ). 

The possible existence of IMBHs in globular 
clusters is suggested by the Mbh — cr relation 
for galaxies, where Mbh is the mass of the cen- 
tral massive BH, a the velocity dispersion of the 
bulge. Extending this relation down to veloc- 
ity dispersions typical for globular clusters (« 
lOkms^^), we expect BH masses in the range 
of ^ 10"^ — 10'' Mq. There are several claims 
in the literature that, indeed, such IMBHs are 
present in some globular clusters. The evidence 
is mainly based on measurements of the stellar 
velocity distribution, with velocity dispersions in- 
creasing strongly towards the center if an IMBH 
is present. The most promising candidate s are 
the clus ters Gl (ICebhardt et al.ll2005l 120021 ) and 
w Cen (jNoyola et al .I l20nj^ . However, the pres- 
ence of IMBHs in these clusters is still debated. 
Fo r instance, based on scale d direct TV-body mod- 
els iBaumgardt et al.l (|2Q03f) show that there is no 
need to invoke the presence of an IMBH in the 
center of Gl in order to explain the observed 
velocity dispers i on pr ofile. On the other hand. 



Gebhardt et al.l (|2005l ) point out that, since the 



region where the IMBH influences the velocity dis- 
tribution is barely resolved, additional information 
from higher order velocity moments must be con- 
sidered. Using orbit-based models to fit both sur- 
face brightness and velocity data simultaneously, 
they find, similar to their earlier results, that only 
these higher order velocity moments provide suffi- 
cient evidence for an IMBH in Gl with a mass of 
1.8(±0.5) X 10^ Mq. However, these results must 
still be confirmed using fully self-consistent evolu- 
tionary models that do not rely on assumed mass- 
to-light ratio profiles. 

Numerous ULX sources have been identified by 
Chandra and XMM-Newton, often associated with 
starburst environments. These sources have X-ray 
luminosities of Lx > lO'^^ergs"^, exceeding the 
angle-averaged fiux of a BH of mass < 20 M0 ac- 
creting material at the Eddington luminosity, Le- 
Although many ULX could be identified as active 



galactic nuclei at the centers of galaxies, and thus 
SMBHs, some are clearly off-ce nter with typical 
proje cted distances of « 400 pc (jMiller fc Colbert 
2004f) . For instance, the ULX X41. 4-^60 in 



the s tarburst galaxy M82 (jPortegies Zwart et al 



20041 ) has been found 7" away from the galaxy cen- 
ter which correspon ds to a distance of « 200 pc. 



Kaaret et al.l (|2001l ) argue, therefore, that this X- 



ray source is unlikely to be an under-luminous 
SMBH, as dynamical friction would cause it to 
spiral into the nucleus of the galaxy on timescales 
much shorter than the age of the galaxy. Based 
on this argument they derive an upper mass limit 
of « 10^ M0. From the total X-ray luminosity 
and the assumtion that the ULX so urce radiates 
isotropically at the Eddington rate, iKaaret et al.l 
(I2OOII) also derive a lower limit of 500 M©. How- 
ever, they point out that this value can be lower 
by a factor of a few when mild beaming is consid- 
ered, thus bringing it closer to the stellar mass BH 
range. On the other hand, new er Chandra obser- 



vations by iKaaret et al. ( 2009[ ) have shown that 
during an outburst, where X41. 4-1-60 increased 
its X-ray luminosity by a factor of more than 
3, its X-ray spectrum remained in the so called 
"hard state" , characterized by a dominant power- 
law component containing 80% of the total flux 
(|Remillard fc McClintoc^ I2OO6I ). From stellar- 
mass BH X-ray binaries it is known that their 
spectra are in the hard state if i < O.Si^;. If one 
assumes the same limit for higher-mass BH X-ray 
binaries, the peak luminosity of 8.5 x 10^°ergs~^ 
inferred from Chandra observations would imply 
a lower mass limit of « 2000 Mq (|Kaaret et al 



20091 ). As one can see, even if one considers that 



the radiation could be mildly beamed, the mini- 
mum mass is well above the stellar-mass BH range, 
making this source a prime IMBH candidate. An- 
other good candidate for an IMBH is the recently 
found ULX source HLX-1 in the edge -on spiral 
galaxy ESO 243-49 (|Farrell et al.ll2009l ). With a 



maximum luminosity of up to 1.1 x 10 ergs in 
the 0.2 -10 key band, and a ccounting for beaming 
effects, iFarrell et all (|2009f) obtain a conservative 
lower limit of ~ 500 Mq. 

Further indications for the presence of IMBHs 
in some globular clusters come from unusually 
high mass-to-light ratios measured in their centers, 
as inferred from pulsar timing measurements. For 
instance, the Galactic globular cluster NGC 6752 
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has five millisecond pulsars; t hree of which are in 
the core ( D'Amico et al ][2Qq2). Two of these three 
have negative period derivatives and one has an 
anomalously high positive period derivative. If the 
spin derivatives are due to the gravitatio nal poten- 
tial of the cluster, iFerraro et al. (I2OO3I) conclude 
that the mass-to-light ratio in the core is M/L 
PS 6 — 7, much higher than inferred for most glob- 
ular clusters. For comparison, M/L « 2 is what 
one would expect for an old star cluster wi th stan- 
dard IMF based on stellar evolution alone (jCaputo 
19851 ). Furthe rmore, based on the position of one 
of the pulsars. iFerraro et al.l (|2003l ) find that there 
is 1000 — 2000 Mq of underluminous matter within 
the inner 0.08 pc of the cluster. This could be ex- 
plained by a ~ 1000 Mq IMBH in the center of the 
cluster, but also by an exceptional concentration 
of dark remnants, or a ~ 100 Mq black hole that is 
offset but near the projected location of the three 
millisecond pulsars in the cluster core. As the high 
spatial resolution in their o bservations does not 
show a cusp down to 0.08 pc, Ferraro et al. ( 20031 ) 
conclude that any central BH must have a mass 
M < 1000 Mq. 

There has been s ome c onsiderable theoretical 
work by I Gill et aL ( 2008 ) on observational evi- 
dences of IMBHs in star clusters. Mass segre- 
gation in star clusters, brings heavier stars to- 
wards the centre increasing the average stellar 
mass in the core. As a diagnosti c tool to quan- 
tify mass segregation. Gill et al. ( 20081) defined 
this radial variation of the average stellar mass as 
Am = {m}^ — where Am is the difference 

of the average mass in the core ((m)^) and the 
average mass at the projected half mass radius of 
the cluster {{m)^^^). A cluster with no IMBH, on a 
relaxation timescale settles to a quasi-equilibrium 
configuration with varying degrees of mass segre- 
gation (Am). They found that in simulations of 
clusters with an IMBH, mass-segregation (Am) is 
significantly quenched. The idea is that, if there is 
an IMBH in the core, since the IMBH will be more 
massive than any of the other massive stars in the 
core, the IMBH has an extremely high probabil- 
ity of exchanging into a binary in a close three- 
body encounter. The subsequent interactions of 
this IMBH in a binary, with other massive stars in 
the core might kick out or scatter the other mas- 
sive stars in the core. Since it will be the mass- 
segregated massive stars which will be scattered 



out of the core in this way, (m)^ will decrease. 
Thus, in a cluster with no central IMBH, mass 
segregation (Am) will be more pronounced than a 
cluster with an IMBH. According to the authors 
this phenomenon of quenching of mass segregation 
in presence of an IMBH can be observed by high- 
resolution imaging of the of the cores of a large 
sampl e of globular clusters b y Hubble Space Tele- 
scope ( Pasquato et al.ll2009l ). 



1.2. Pathways to IMBH formation 

There are several pathways discussed in the lit- 
erature through which IMBHs may form. The 
simplest way is through the core collapse of a 
massive Population HI star for med in a mini 
dark matter halo at high redshift (jMadau &: Rees 
20011 ). At lower redshifts such massive stars must 
be grown through mergers of lower mass stars, 
which requires ra ther large stellar densities, typi- 
cally ^ lO^pc"^ dFreitag et al.ll2006al: lA"rdi et al 
2008t iBaumgardt et al.l bOOsT As an alterna- 



tive to stellar mergers it is also possible to in- 
crease the mass of a stellar mass BH by tidally 
disrupting and accreting other stars; however, 
this seems to require ev en larger stellar densities 



( Baumgardt et al.l 120081 ) 



In a cluster with a broad range of stellar masses, 
large stellar densities can be achieved through 
mass segregation, which will drive the most mas- 
sive objects to the center, while lighter stars 
spread out to attain kinetic energy equipartition. 
However, for any reasonable IMF the most massive 
stars are unable to achieve energy equipartition 
with the lighter stars and therefore decouple from 
the rest of the cluster and form a compact subsys- 
tem in the cluster center. This process is called 
Spitzer instability ( Spitzer 19691) and causes core 



collapse to be accelerated (lSpitzeiill96S:lVishniac 



1978i lWatters et al.lliooi ICiirkan et al 
if (I2OO4L 



shown by iGtirkan et al.l 



2004h . 



As 

hereafter Paper 
I), for a realistically broad stellar mass spectrum 
the core collapse time (tec) can be as short as 
tec = 0.15trc(0), where trc(O) is the initial central 
relaxation time. 

Once the most massive objects have segregated, 
the formation of an IMBH can occur in two dif- 
ferent ways: one is when the massive stars have 
already evolved into stellar-mass BHs at the time 
of core collapse and these BHs then merge by 
emitting gravitational waves. The BHs first form 
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BH binaries which then harden through dynami- 
cal interactions with other objects until the binary 
is close enough for gravitational radiation to dissi- 
pate sufBc ient orbital energy u ntil the BHs merge 
(see, e.g., O'Learv et al. 2006[ ). However, grow- 
ing an IMBH with a mass of ~ 1000 M0 through 
such stellar-mass BH mergers is only realistic for 
the very massive clusters such as those in galac- 
tic nuclei ( O'Learv et al. 2006[l. For smaller sys- 
tems s uch as globular clusters iMouri fc Taniguchi 
(I2OO2I ) and lO'Learv et all (|2006l ) have shown that 
this mechanism is rather inefficient. This is be- 
cause it is much more likely that the stellar mass 
BHs escape through strong dynamical binar; 



200' 



interactions (|Portegies Zwart fc McMillan 
or recoil from asymmetric gravitational wave 
emiss i on dMiller fc Hamilton 2002 ; Giiltekin et al 



20041 : lO'Learv et al.ll2007h given the low cluster 
escape speed. 

Another way to form an IMBH, which is not 
restricted to galactic nuclei, is through mergers 
of massive main-sequence stars that segregate to 
the center and drive cluster core collapse before 
formation of stellar BHs. This subsystem of mas- 
sive stars can enter a phase of rapid collisions, 
and since the most massive object has the largest 
cross-section for further collisions, this object is 
expected to grow in a runaway fashion. The re- 
sulting V MS event ually collapses to for m an IMBH 
jPortegi es Zwart & : McMillan 2002; Gii rkan et al. 
|2004| ). As the time it takes for the most mas- 
sive st ars to turn in to BHs is approximately 
3Myr (|Mevnet fc Mae dcr 2000), an IMBH is only 
formed through runaway merging if the cluster 
reaches core collapse within the first 3 Myr of dy- 
namical evolution. 

This simple picture has a few important 
caveats. First, the fate of such a massive merger 
remnant formed by a runaway is rather uncertain. 
Direct monolithic collapse to a BH with no or lit- 
tle mass loss is a possible out come, at least for 
suffic iently small metallicities (jHeger fc Wooslev 
l2002l) . However, it has been suggested that mass 
loss from stellar winds could dominate the mass 



increase due to repeated mergers (.Glebbeek et al 



20091 ). In this case it might be difficult to form a 
VMS at all. On the other hand, within the run- 
away phase, it has been shown that, for clusters 
like those studied here, the time between colli- 
sions is much shorter than the Kelvin-Hclmholtz 



timescale of the collision product, so that the 
growing VMS must be out of ther mal equilib- 
rium. Instead the calculations of iGlebbeek et al 



( 2009t ) were done assuming that each merger rem- 
nant behaves exactly like an ordinary massive star 
in thermal equilibrium. The uncertainty in the fi- 
nal VMS (and IMBH) mass is further increased 
when considering that a stronger wind mass loss 
would also lead to a stronger expansion of the 
cluster core, decreasing the mass growth rate of 
the VMS by lowering the collision rate. To ad- 
dress this problem, one would need to perform a 
fully self-consistent simulation coupling the stellar 
dynamics with detail radiation hydrodynamics of 
the stellar collisions and mass loss from merger 
remnants. This is clearly beyond the scope of this 
paper. 

The runaway collision scenario in globular clus- 
ters h as been extens i vely in vestigated in Paper I 
and in iFreitag et al. I (l2006allbh using Monte Carlo 
(MC) simulations for a large variety of initial con- 
ditions. Paper I focuses on the dependence of t^c 
on the shape and the width of the IMF, the pres- 
ence of a Galactic tidal field, and the cluster den- 
sity profile. The key result is that for clusters with 
a broad range of masses, ice is set by the cen- 
tral relaxation time, trc(O), which means that for 
multi-mass clusters core collapse depends on the 
local conditions in the core while for single-mass 
clusters core collapse is a global phenomenon. Fur- 
thermore it is clear that ice depends only weakly 
on the external tidal field. It was also found that 
the dependence of tec on the mass spectrum can 
be conveniently expressed by a single parameter 
rrimax/ {m), where iTT-max is the maximum and {tti) 
the average stellar mass. For rrimax/ ("i) > 50 the 
ratio icc/iic converges to a constant value w 0.15 
for all IMFs and cluster density profiles. Com- 
bined with the requirement that tec be less than 
3 Myr, this relation provides a uniform criterion 
for runaway growth to occur for a large range of 
possible, unsegregated cluster configurations. 



Freitag et al.l (|2006al ) performed similar simu- 



lations but also incorporating collisions explicitly. 
They quantified the dependence of the onset of the 
runaway on the initial collision time, and found 
that runaway growth happens earlier with respect 
to irh, the half-mass relaxation time for initially 
more collisional clusters. Thus, collisions extend 
the parameter space of initial cluster conditions 
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for runaway to occur. However, as pointed out 
by the authors, for any standard IMF (Kroupa, 
Salpeter, Miller-Scalo) this effect is neghgible for 
masses typical for globular clusters. Thus, in this 
regime, the condition for runaway to occur reduces 
to the one found in Paper I, based on the central 
relaxation time alone. 

The runaway collision scenario has also been 
ver ified numerically by direct iV-body si mulations . 
In IPortegies Zwart fc McMillanI (|2002l) runaway 
collisions were produced in sufficiently dense and 
highly concen trated clusters with on l y a fe w 10'' 
stars initially. IPorte gies Zwart et al.l (|2004 ) mod- 
elled the evolution of MGG-11 with w 10^ and, 
in one case, with « 5 x 10^ stars, and found 
that, similar to Paper I, only clusters with a 
short enough mass-segregation timescale, or, cor- 
respondingly, low enough irh, are likely to pro- 
duce a runaway object. However, their results also 
show that, in addition to a short irh, these clus- 
ters must also be sufficiently concentrated, cor- 
responding to King models with Wq > 9, in 
order to trigger a runaway object. This might 
be due to the fact that in the iV-body simula- 
tions binaries formed by three-body interactions 
(jFreitag et al.ll2006a[) . a process which is not in- 
cluded in the MC runs, and these binaries dom- 
inated the coUisional evolu tion in the core. On 
the other hand, for TV > lO'^ lFreitag et al.l (|2006ah 



have demonstrated that three-body binaries are 
of littl e importance for the c o llision process. Sim- 
ilarly, IPortegies Zwart et al. I too4 find for their 
large N ^ 5 x 10^ run, that the influence of three- 
body binaries on the collisional evolution is lower 
than in their lower- iV runs. Thus, it appears that 
the nature of collisional runaway growth in a star 
cluster changes with increasing N from being less 
dominated by binary collisions and involving more 
and more single-single collisions. 

Although, the runaway merger scenario has 
been extensively investigated for a large variety 
of initial conditions, almost all of these studies 
started with unsegregated clusters, that is, clus- 
ters, where the stellar mass is not correlated with 
the radial position within the cluster. However, 
in recent years observations frequently indicated 
that many young star clusters with ages of only 
a few Myr show already a strong degree of mass 
segregation, suggesting that mass segregation may 
be primordial and motivating a new study of its 



effects on collisional runaways. 

1.3. Primordial Mass Segregation and 
Top-Heavy IMF 

Significant mass segregation has been found in 
many young star clusters. Mass segregation can be 
understood as the tendency towards equipartition 
of kinetic energies. This tendency for equipartition 
is a consequence of gravitational encounters, which 
attempt to drive the local velocity distribution 
toward a Maxwellian , with mi (Vi) = m2 {V2) 
(jHeggie fc HutI l2003l Chapter 16). As a result, 
massive stars move more slowly, on average, than 
lighter ones, so the massive stars drop lower in 
the potential well, while the stars o f smaller mass 
move out and may even escape ( Heggie fc: Hull 
20031 . Chapter 16). This results in having the more 
massive members of a gravitationally bound sys- 
tem closer to the center whereas the lighter mem- 
bers are found further away. This dynamical mass 
segregation acts on a time-sca le of the order of the 
relaxation tim e of the cluster ( Heggie fc Hutll2003l: 
Spitzeij[T987ai Section 4.2). 

However, there are indications that the degree 
of mass segregation seen in many young clus- 
ters cannot be a result of their dynamical evo- 
lution as their ages are much less than their relax 
ation time (Hil lenbrand 1997; Fischer et al. 19£ 
de Griid Il998l: IffiUenbrand fc HartmannI 1 199 



Gouhermis et al.ll2004 IStolte et al.ll2006l ). For in- 
stance, for NGC330, the richest young star cluster 
in the small Magellanic Cloud (SMC), observa- 
tions show that the IMF becomes steeper at in- 
creasing distances from the cluster center, with the 
number of massive stars decreasing from the core 
to the outskirts of the cluster 5 times more rapidly 
than the less-massive objects, while the age of 
NGC330 is 10 times shorter than the expected re - 
laxation time of the cluster (jSirianni et al.l 120021 ). 
Another example is provided by the Orion nebula 
Cluster, in which it can be argued from numer- 
ical results that the massive stars in the centre 
cannot have formed in the outer regions. This 
implies that the stellar mass is to some degree 
a fu nction of the initial posi tion within the clus- 
ter ( Bonnell fc Davie^ 19981) . Interestingly, the 
3 mm continuum observations obtained with the 
Combined Array for Research in Millimeter- Wave 
Astronomy (CARMA) of 11 Infrared Dark Cloud 
(IRDC) cores establishes that mass segregation 
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can be ide ntified at the formative stage of a stel- 
lar cluster (jPerez Munoz fc Carpenteiil2007l ). All 
these observations suggest that many clusters may 
have been formed significantly mass segregated. 

Primordial mass segregation has been explained 
either by the formation of massive stars preferen- 
tially in t he densest regions of the parent molecu- 
lar cloud ( Murray fc Lin 1996') or by competitive 
gas accr etion during the earlies t phases of star for- 
mation ( Bonnell fc Batdl2006l) . An additional al- 
ternative has been suggested by i McMillan et al 



(|2007l ). They show that a much higher degree of 
mass segregation in a cluster can be achieved if the 
cluster is a result of one or more mergers of smaller 
clusters. As the smaller clusters have shorter re- 
laxation times, mass segregation proceeds faster. 
When the clusters merge, this larger degree of 
mass segregation is preserved in the final cluster. 

Primordial mass segregation also has important 
consequences. It has been shown that, in an ini- 
tially mass-segregated cluster, the effect of early 
mass loss due to stellar evolution is, in general, 
more destructive than for an unsegregated cluster 
with the same density profile, since this leads to 
shorter lifetimes, a faster initial evolution toward 
less concentrated structure and a flattening of the 



stellar mass function ([McMillan et al.ll2008f) . 

Recently, lArdi et all (|2008l ) studied the influ- 
ence of initial mass segregation on the runaway 
growth of a massive star by means of direct iV- 
body simulations of up to ~ 131 stars . Cont rary to 
the expectations from iGiirkan et al. 1 t004\ . they 
found that, for a given density proflle, initial mass 
segregation does not increase the available param- 
eter space of cluster initial conditions leading to 
runaway growth. They argue that this is be- 
cause initial mass segregation decreases the col- 
lision rate of stars in the core as, due to the in- 
creased average mass, the number density is de- 
creased. This is in line with earlier TV-body simu- 
lations suggesting that runaway growth can occur 
only wh en the cluster is initially su fficiently col- 
lisional, (|Portegies Zwart et al.l 120041 ) . in contrast 
to predictions based only on the core collapse time. 
This result should be tested for larger N as the na- 
ture of the collisional runaway changes from being 
dominated by three-body binary formation at low 

to single-single collisions at high A^. 

One of the main uncertainties in star clus- 
ter evolution lies in determining the true initial 



mass function (IMF). Often it is assumed that 
the IMF is a standard power-law (or power-laws 
with different indices in different mass ranges) 
with no primordial radial variation in the clus- 
ter (e.g. . [Salpeteil [1959I: iMiller k Scalol Il979a : 



Kroupal 120011) . Deviations from the standard 
IMFs are observed in many cl usters both at th e 
high and low- mass ends (e.g., Elmegreen 2004 ). 
In particular at the high mass end MFs are ob- 
served to be generally flatter compared to stan- 
dard Salpeter power-law in y oung massive clus- 



ters like the A rches cluster (jStolte et al, 
Kim et aLlliooj ). 



2002; 



Numerical Methods 



2.1. Monte Carlo code 



The numerical method that has been used here 
to investigate the dynamical evolution of star clus- 
ters is the MonteCarlomethod, based on the clas- 
sic work_of Hino3( 1971 ) and described in detail 



Fregeau fc Rasiol ( 2007 , and references therein) . 



In Monte Carlo simulations, A^, the total number 
of stars in the cluster is dependent on the initial 
half mass relaxation time in t he cluster. Fo r a 
Plummer sphere it is given by ( Spitzer 1987bl eq. 
2.63), 



irh(O) 



0.138Af 



In-fcN V GM J 



1/2 



(1) 



where 7c ^ 0.01 is the coulomb logarithm and 
is the half mass radius. 

Since our code now includes an explicit treat- 
ment of all stellar collisions we briefly summarize 
here t he 'sticky sphere' me thod for stellar colli- 
sions (|Freitag fc Benzll2002t ). 

In the sticky sphere approximation a collision 
occurs whenever the centers of two stars pass 
within a distance d = (i?i + i?2), with Ri^2 being 
the stellar radii. Until this distance is reached, 
the gravitational influence of other stars as well 
as any mutual tidal interactions are neglected. 
The cros s section for such a collision is given by 
Binnev fc Tremaine (1987. section 7.5.8), 



S. 



12 



12\2 



^(i?i+i?2)'(l + %^), (2) 

Kel 



where bmax is the largest impact parameter lead- 
ing to contact, Vrei is the relative velocity between 
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two stars and vl"^ = {2G{Mi+ M^) / {Ri+ R2)f -^ ■ 
In a cluster where all stars have the same mass 
M* and radius Rt, , the average local collisi on time 
Tcoii is given bv lBinnev fc Tremaind ( 1987 , section 
7.5.8) 



1 — / d?VidPv2f{vi)f{v2)\\v^ - V2\\S, 



(3) 

where is the stellar number density. For a 
Maxwellian velocity distribution this becomes 



T,, 



= 16\/T:nf.ayRl 



1 



2a?, R^ 



(4) 



where (T„ is the velocity dispersion. 

In order to resolve collisional processes in a MC 
simulation we constrain the time step size St ac- 
cording to an estimate of the central collision time 
using 

St < ff,M, (5) 

where / = 5 x 10"'^ is a constant chosen small 
enough to ensure that collisions are sampled suf- 
ficiently, and TcoU is an estimate of Tcoii based on 
equation SI and given by 



1 



coll 



1 



GjAUR^) 



(6) 



with quantities in angular brackets being local av- 
erages. The collision probability for two neighbor- 
ing stars is calculated as 



Pr, 



n^VrelScollSt, 



(7) 



where n* is a local estimate of the stellar number 
density. The time step size St, is chosen such that 
Pcoii << 1. A random number (between and 
1) is then drawn. If the random number < Pcoii 
two stars are selected for a collision and they are 
merged under the assumption of mass and linear 
momentum conservation. 

The units adopted for our simulations are the 
standard iV-body units of the MC code, and the 
same as in Paper I. 

2.2. Implementation Of Primordial Mass 
Segregation and Top-Heavy IMFs 

The effect of primordial mass segregation was 
studied in paper I using a preliminary prescrip- 
tion. However, the recipe in Paper I is not defined 



with well constrained parameters. In this paper 
we have used recipes which are more realistically 
modelled and we do a much more extensive study. 
Here we consider with two prescr iptions for 



mass segregation: (i) the prescription of lSubr et al 



(|2008i) for generating mass-segregated clusters, 
in which the degree of mass segregation can be 
adjusted by a parameter S related to the mean 
inter-particle ene r gy of the stars; and (ii) the 
iBaumgardt et al" 1 (l2008l) prescription, which cre- 
ates a maximally mass-segregated cluster in virial 
equilibrium. 

From TV-body models ISubr et all (|2008[) find 
that the quantity that is transferred between the 
light and massive stars is the potential energy, 
while their average kinetic energy remains nearly 
constant during the cluster evolution. Therefore, 
mass segregation is generated in terms of mean 
inter-particle potentials (C/*-' ), which is parame- 
terized as 



ML 



ML 



(8) 

where and M^j denote the masses of the ith 
and jth particle, Mtot denotes the total mass of 
the cluster, C/tot is the total potential energy of the 
cluster, Mgj^^jj is the sum of all M^k < M^ij and S 
is the degree of mass segregation which can have 
values between and 1 (for a de tailed description 
of the code refer to ISubr et all ([2008)). In this 
recipe S* = implies an unsegregated cluster and 
5 = 1 a completely mass segregated cluster. Re- 
lating to entropy S — has the lowest entropy 
and the system is highly symmetric in terms of 
< >, while for S = 1, all the binding en- 
ergy is carried by the two most massive stars in 
the cluster and corresponds to a state of maxi- 
mum entropy. With this parameterization of the 
inter-particle potential, quasi-stationary, star-by- 
star realizations of mass-segregated clusters are 
generated. The stars are assigned masses accord- 
ing to the Salpeter IMF with maximum and min- 
imum mass of 0.2 M0 and 120 M© respectively. 

An important property of this mass segregation 
recipe is that both the density profile and the av- 
erage mass in the core ((m)^) is different for each 
value of 5*, starting from a Plummer sphere for 
5" = 0. This has the consequence that irc(O) which 
is dependent on density in the core (pc) and (m)^ 
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changes. This in turn causes tec to change. 



(|2008h 



on 



The recipe from iBaumgardt et al 
the other hand, does not change the underlying 
density profile and generates mass-segregated clus- 
ters, in comparison, much faster. It essentially 
sorts all stars such that, for a given density pro- 
file, the most massive stars have, on average, the 
lowest specific total energy. 

Figures [T] and [5] show the average mass profile 
for the initially mass-segregated clusters (and un- 
segregated c luster s also for comparison) us i ng the 
Subr et al.l (|2008l) and IBaumgardt et all (|2008l ) 
recipes, respectively. As can be seen in these 
plots, the average mass rises rapidly within one 
half-mass radius, by up to a factor of ~ 5 (for 
the Baumgardt et al. (2008) recipe as well as for 
S = 0.3) for the initially mass-segregated clus- 
ters. Similarly Figures [3] and H] show the enclosed 
mass for the same two prescriptions. Compar- 
ing the median positions of the massive stars in 
these models, it becomes clear that the massive 
stars are on average at much smaller radii in ini- 
tially mass-segregated clusters compared to non- 
segregated ones, as expected. The median dis- 
tance from the cluster center can be up to a factor 
~ 4 smaller for stars with mass larger than 50 M0 . 
Comparing the enclosed stellar mass profiles for 
the two recipes ( Figure El and Figure Bt it can 
be seen that, the IBaumgardt et al.l (l2008f) recipe 
agree pretty closely with the Subr et all ( 20081 ) 
recipe for S — 0.3. 

In this paper we also study the effect of flatter 
IMFs on the runaway collision scenario in young 
massive star clusters. We introduce a new IMF 
in this pa p er mo tivated by the Arches IMF from 



Dib et al.l (|2007[ ). For the rest of the paper, this 
new IMF will caUed the "Variable IMF" and we 
will denote the sections of this IMF {dN/dm (x 
m~") as [ai . . .ai] with ^4 corresponding to the 
high mass tail of the IMF. Th is Variable IMF is a 
variation on the Arches IMF (jPib et al.ll2007t ) by 



leaving out the middle section from 3.0 — 15.0 Mq 
so that the IMF is then more easily compared 
to the standard Kroupa IMF. The upper mass 
section of the Variable IMF is much flatter than 



traditional ISalpeteij ( 19551 ) . iKroupal (|2002l ). or 
Miller fc Scald (|l979bl ) IMFs. We therefore ex- 
pect a greater number of high mass stars, increas- 
ing both the average mass as well as the entire 
cluster mass. 
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Fig. 1. — Average stellar mass profile for initially 
mass-segregate d clusters created according to the 
prescription of Subr et al. ( 20081 ). All the simu- 
lations shown here contain 5 x 10^ stars with a 
Salpeter IMF (within Mi„i„ = 0.2 Mq and M^ax = 
120 M0). The horizontal axis shows the cluster ra- 
dius in units of the initial half-mass radius. The 
dotted straight line is for an initially unsegregated 
cluster and the red and black solid straight lines 
are for clusters with different degrees of primordial 
mass segregation (5 = 0.3, 0.2). The average mass 
for the mass-segregated clusters increases steeply 
within ^ 1 rii. For S — 0.3, the increase is by up 
to a factor of ^ 5. 



2.3. Initial Conditions 

For all the simulations in this paper our code in- 
cludes physical processes such as two-body relax- 
ation and physical stellar collisions between stars 
in the 'sticky sphere' approximation. In this paper 
we do not include stellar evolution in our calcu- 
lations since our aim here to investigate dynam- 
ical processes taking place before even the most 
massive main sequence stars in the cluster have 
evolved (i.e., < 3Myr). We have also not included 
any primordial binaries in our simulations just to 
keep the simple picture of runaway in presence of 
primordial mass segregation, and exclude the ef- 
fects of primordial b inaries on runaway collisions 



(iGiirkan et al.l|2006l ) 



The most important physical properties of all 
our initial cluster models are given in Table [2] and 
Table [3l A comparison of the variable IMF with 
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Fig. 2. — Same as Fig[T] but the primordially 
mass-se gregated cl uster is generated using the 
feaumg ardt et aL (I2OO8I) recipe. Even in this 
case the average stellar mass of the cluster rises 
strongly within 1 Th, by up to a factor of '-^ 5. 



the Arches IMF or with standard Kroupa IMF is 
given in Table [TJ In Table [2] we have listed simula- 
tions of clusters with the Variable IMF, varying 
and a4. We have also listed the simulations done 
with Salpeter IMF varying Mmax- In all these sim- 
ulations our focus has been to investigate the effect 
of IMF on core-collapse of young massive clusters. 
All our models start with an isolated Plummer 
sphere and the virial radius has been chosen such 
that the corresponding cluster has a core-collapse 
time of 3 Myr. 

In Table [3] we have listed simulations of clusters 
for different values of initial N and initial virial 
radius, r-vir, varying the degree of primordial mass- 
segregation. Virial radius defined by 



GM2 



(9) 



where M is the total mass and E is the total grav- 
itational energy of the cluster. This grid formed 
by the different values of initial N and initial rvir 
is referred to as the parameter space in the later 
part of the paper. For all the simulations listed 
here our aim has been to investigate the effect of 
primordial mass-segregation on ice of young clus- 
ters with different initial conditions {N and rvir). 
We have calculated ice based on the central relax- 



ation time ire defined as 



trc(O) 



^7rG2(ln7c7V)n(TO) 



(10) 



where CTc, n, {m)^ are the three-dimensional 
velocity dispersion, number density, and aver- 
age stellar mas s respectively at the cluster center 
( Spitzeij Il987al eq. 3.37). As a typical refer- 
ence model similar to Paper I, all the simulations 
in Table [3] are initiated with Plummer models 
and Salpeter IMF (within Mmin = 0.2 Mq and 
Mmax = 120 M0). A nother reason for u sing Plum- 
mer sphere is that the lSubr et al. ( 2008 ) formalism 
starts with an isolated Plummer sphere. 

3. Results 

3.1. Primordial Mass segregation 



Freitag et al.l (|2006af ) were the first to study the 
core collapse and collisional runaway for unseg- 
regated clusters using a MC simul ation code in- 



cluding stellar collisions. Similar to iFreitag et al 



(j2006al ) we did not take into account stellar evo- 
lution in the simulations since all the simulations 
were limited to the first 3 Myr, before the most 
massive stars lose mass in supernovae explosions. 
If a cluster had a core collapse time more than 
3 Myr we implicitly took stellar evolution into ac- 
count by ending the simulation at 3 Myr. Since 
our code now includes stellar collisions, we first 
checked whether we are able to reproduce their 
results. They found that for all models with 
trc < 20 Myr, runaway formation of a VMS oc- 
cured. 

Figure [51 is a pa rameter survey similar to 
Freitag et al.l (|2006al their Fig. 1), showing 
for each simulation whether a runaway occurred 
(filled circles) or not (open circles), for all our 
unsegregated models varying r^ir and the num- 
ber of stars in the cluster. The solid straight line 
corresponds to tec = 0.15trc — 3 Myr, not includ- 
ing collisions. Simulations of clusters with initial 
conditions lying below the straight line will have 
tec < 3 Myr, whereas tec > 3 Myr for simulations 
with initial conditions above the line. It can be 
clearly seen in Figure [5] that the initial conditions 
leading to a runaway fall indeed below this line. 
T hus we have successfu lly reproduced the results 
of lFreitag et"aD (|2006al ) with our code and have 
reconfirmed the validity of the simple criterion for 
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runaway collisions to occur in young dense star 
clusters. 

We then repeated the same set of simulations 
for primord ially mass-segrega ted clusters using the 
recipe from lSubr et al.l (|2008l ) with S > 0. Figure 
[6]shows our results for mass-segregation parameter 
5' = 0.30. We see that, with primordial mass seg- 
regation, the initial cluster r^ir can be chosen sev- 
eral times larger compared to unsegregated clus- 
ters and still lead to a runaway. Thus, simulations 
of primordially mass-segregated clusters increase 
the parameter space (range of rvir in pc) for run- 
away collisions to happen when compared to sim- 
ulations of primordially unsegregated clusters. 

To further illustrate the effect primordial mass 
segregation has on the mass growth of the most 
massive star we compared results of a primor- 
dially mass-segregated vs an unsegregated cluster, 
with otherwise identical initial conditions. Fig- 
ure [7] shows the growth curve of the most mas- 
sive star in an unsegregated cluster and in a pri- 
mordially mass-segregated cluster {S = 0.3), with 
iV = 5 X lO'^ stars and Tvir — 0.64 pc. The unseg- 
regated cluster does not have a runaway or a very 
steep mass growth before it reaches core collapse, 
and only a 200 Mq star is formed within 3 Myr. 
For the primordially mass-segregated cluster, we 
clearly see a very steep mass growth leading to 
a formation of a 900 Mq star within 3 Myr and 
henc e a runaway. A gene ral result, in agreement 
with iFreitag et al.l (|2006al) . is that only one VMS 
forms in the cluster and there is no sign of multiple 
runaways. Note that this result might be different 
for clust ers with significant f ractions of primordial 
binaries (jGiirkan et al.ll2006[ ). As in this example, 
we also note in all our simulations that the time 
of the runaway coincides with the time of core- 
collapse, which decreases from t^c ^ 6 Myr for the 
unsegregated cluster, and to tec ~ 2 Myr for the 
primordially mass-segregated cluster. 

To show how strongly the core collapse time de- 
creases for other values of S we plot in Figure[8]the 
core collapse time for clusters with different initial 
Tvir, against S. We clearly see a trend of decreas- 
ing tec with increasing S. The decrease in tec is 
steepest for 5 < 0.1, causing a reduction of tec by 
a factor of « 2. For S > 0.1 the decrease in tec is 
weaker, reducing tec by another 30%. This trend 
is also illustrated in Table |3] where we have shown 
that tec decreases for simulations with primordial 



mass-segregation, when compared to simulations 
of initially unsegregated clusters with similar ini- 
tial conditions. 

To analyze this trend further we first check 
whether the simple relation from Paper I between 
tec and tic(O) (tec/trc(0) — 0.15) remains still valid 
for mass segregated clusters. In Figure |9] we plot 
tcc/trc against S for clusters with different Tvir. We 
see that tcc/trc remains nearly constant at 0.15 
with only a ^ 10% scatter. Thus, we conclude, 
that the ratio of tcc/trc remains consistent with 
the value found in Paper I and the strongly de- 
creasing tec can then only be caused by a decrease 
in trc(O). This is also shown in Figure [101 where 
we plot tie(O) against S for clusters with different 
Tvir. We notice a very similar decrease in tj-c as in 
tec, with increasing S, so the reason for the shorter 
tec must be attributed to shorter tic. 

However, in the Subr et al. ( 20081 ) formulation 
it is not a priori clear whether this decrease in t^c is 
entirely related to an increase in primordial mass 
segregation, as the central density (pc) also in- 
creases with S (§4.2). In order to disentangle these 
effects, we calculate relative contributions of all 
the factors to tj-c. We note that t^e oc a^/pc {m)^ 
where ac is the velocity dispersion in the core, and 
(m)^ is the average stellar mass in the core. Fig- 
ure [Tl] plots the variation of a^, 1/pc and 1/ {m)^, 
normalized to their values at S" = 0, against S in 
a cluster with = 5 x 10^ stars. The inverse of 
the normalized initial central density of the clus- 
ters decreases as S is increased from 0.0 to 0.15 by 
25%, and then remains nearly constant for larger 
S values. The inverse of normalized (m)^ in the 
cluster always decreases with the increase in mass 
segregation, up to a factor of 2.5 for S — 0.3. The 
variation of crj?, on the other hand, is only by a 
factor of w 10%, implying that Cc stays almost 
constant in the core and does not change with S. 
From this it follows that for S < 0.15 both the in- 
crease in (m)^ as well as pc contributes equally to 
the decrease in t^c. On the other hand, for larger 
S, tfc is dominated by {m) ^ while the contribution 
from Pc is minimal (w 20%). So it is the increase 
in {m)^ that mainly causes the low tic values for 
larger S. Thus the maximum increase in the pa- 
rameter space for a runaway to occur is driven 
mainly by "true" primordial mass segregation. 

The advantage of including stellar collisions in 
our code is that we are able to directly determine 
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how much mass eventually ends up in the VMS. 
Paper I showed that, as the core collapse proceeds, 
the mass contained in the collapsing core which 
forms the mass reservoir for the runaway converges 
to a value M^c ~ 0.001 - 0.002 Mtot, where Mtot 
denotes the total mass of the cluster. In general 
we find that the fraction of Mtot that ends up as 
a VMS is in a similar mass range as Mcc in Paper 
I. In Figure [12] we show that fraction as a function 
of S for a cluster with = 5 x 10^ stars. It can 
be clearly seen that the fraction of Mtot ending 
up as a VMS increases almost by a factor of 3 for 
S = 0.3 and hence mass of the VMS increases sig- 
nificantly with primordial mass segregation. Thus 
primordial mass segregation not only increases the 
parameter space for runaways to occur, but also 
produces more massive VMS. 

3.2. Top-Heavy IMFs 

As discussed previously, there are indications 
that some young massive star clusters may be born 
with IMFs flatter than the Salpeter IMF. Here we 
study how such flatter IMFs influence the avail- 
able parameter space of cluster initial conditions 
for a runaway to occur. One crucial result from 
Paper I is that the ratio of maximum to aver- 
age stellar mass in the IMF is the most important 
parameter setting the timescale for the onset of 
core collapse. From Paper I if mmax/ {mn) < 40, 
tec oc (mmax/ (to))"^'^ and if mmax/ (m) ~ 50, 
a domain is reached by any realistic IMF where 
tec — O.lS^rc- A flatter IMF increases the num- 
ber of massive stars, thus increasing the average 
stellar mass in the cluster, which in turn causes 
an increase in tec- This means that, although for 
flatter IMFs we have more massive stars, the pa- 
rameter space for runaway collisions to occur ac- 
tually gets smaller. On the other hand, collisions 
can r educe trr. by constantl y dissipating orbital en- 



ergy ( Freitag et al.|[2006al) . This means that the 
increase in tec seen in simulations without an ex- 
plicit treatment of stellar collisions, is to some de- 
gree compensated by the effect of collisions. In 
Figure [13] we quantify this effect for flatter IMFs 
considering a wide range of Wmax/ (w) values. We 
plot tcc/ti-h vs m„iax/ (™) with and without colli- 
sions. As can be clearly seen, with stellar collisions 
turned off, tec increases; but when collisions are in- 
cluded, ice decreases on average by 15%. Overall 
this causes tcc/^rh for a collisional cluster causes 



the core collapse to stay close to the value for 
the standard Salpeter IMF (0.07) all the way to 
mmax/ (to) w 140. 

We also, studied the combined effects of initial 
mass segregation and a flat IMF on the evolution 
of young massive cluster s . Th is is mainly moti- 
vated bv lEspinoza et aL ( 20091 ) where the global 
IMF of the Arches cluster is determined to be 
flatter than Salpeter, while Arches is also known 
to be strongly mass segregated. Note however 
that the uncertainity on the slope is such that 
it may not be incompatible wi th a Salpeter IMF 
and there are models of Arches ( Harfst et al.ll2010l ) 
that can explain the degree of mass segregation 
with a Salpeter IMF and unse gregated init i al con - 



ditions. On the other hand Harfst et al. ( 2010[ ) 



do not rule out the possibility that Arches could 
be primordially mass-segregated. In Figure [T3l we 
show the results f or primordially rnass-seg regated 
clusters using the Baumgardt et alj ( 20081 ) recipe, 
which, unlike the Su br ct al.. (2008,) recipe, can be 
applied to clusters with arbitrary density profiles. 
This way the decrease in tec is only caused by the 
increase in the average stellar mass at the center 
of the cluster. From Figure [T3] we see that, as in 
the unsegregated case, the effect of a flatter IMF 
increases ice, and for mmax/ (to) = 120 is very 
close to the standard Salpeter value. The change 
for TOmax/ (w) = 160 is only by 20%. Note that 
TOmax/ (to) = 120 corrcspouds to either a Salpeter 
IMF with a high mass cut-off at 80 M©, or to an 
IMF with a slope of —2.1 and a full mass spectrum 
(from 0.1 — 120 M©). This latter slope is precisely 
the one suggested for th e global IMF of Arches 
by lEspinoza etldl (|2009l ). If the global IMF of 
the Arches star cluster is also representative of 
the IMFs for much more massive clusters like the 
Westerlund 1, this would imply that the simple 
relation from Paper I can be used to determine 
whether a young massive cluster will undergo a 
runaway or not, without explicitly accounting for 
a different global IMF or primordial mass segrega- 
tion. 

4. Summary and Discussions 

This work is a continuation of our study of 
the runaway collision scenario in young dense star 
clusters. In this paper our goal was to investigate 
the parameter space for runaway collisions to oc- 
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cur, when the effects of stellar collisions, primor- 
dial mass segregation and a globaUy flatter IMF 
(as indicated by observations of young massive 
star clusters) are accounted for. We considered 
clusters having Plummer density profiles, varying 
the initial virial radius, the number of stars in 
the cluster and the IMF. Prim ordial mass segre - 
gati on was generated using t he Subr et al. ( 20081 ) 
and iBaumgardt et al. ( 20081 ) methods. We natu- 
rally expected that primordial mass segregation in 
a cluster should lead to shorter core-collapse times 
since the massive stars start their lives closer to the 
center of the cluster, where the density and there- 
fore collision rates are highest. Indeed, from our 
simulations we found that the core-collapse time, 
tec decreases with increasing the degree of primor- 
dial mass segregation, increasing the initial virial 
radius for runaway collisions to happen, by at least 
a factor of w 3. 

We find that the simple relation between core- 
collapse and central relaxation times, (ice = 
O.lSirc), discussed in Paper I, still holds even with 
primordial mass segregation. The strong decrease 
in tec hence implies a similar reduction in ire, 
which accelerates the core collapse and enlarges 
the parameter space for runaway colli sions to oc 
cur 



This decrease in t,c, in both the Subr et al 



( 2008 ) and Baumgardt et al. ( 2008[ ) prescriptions, 
is caused by the strong increase in the average 
stellar mass in the core. In the S ubr et al.l ( 2008h 
recipe the decrease in tic is also caused by an 
increase in central density, which, however, does 
not contribute much at higher degrees of mass 
segregation {S > 0.15). 

We find that the fraction of the total clus- 
ter mass that is eventually accumulated onto the 
VMS, the possible progenitor of an IMBH is com- 
parable to the mass of the collapsing core, Mcc, 
already determined in Paper I. If IMBHs were in- 
deed formed in massive star clusters, our results 
show that the ratio of the mass of the IMBH to 
the total stellar cluster mass follow the Af/Afbuigo 
ratio. However, the runaway collision scenario im- 
poses the requirement that the IMBH must have 
formed in a cluster with initial central relaxation 
time shorter than 20 Myr. 

Finally we find that for top-heavy IMFs (in un- 
segregated clusters) the parameter space for run- 
away collisions is reduced, since flatter IMFs in- 
crease tec, by a factor w 1.4 for mmax/ {m) w 120. 



However this increase is to some degree balanced 
by collisions, which happen more frequently in 
clusters with flatter IMFs as we get more mas- 
sive stars. In addition, primordial mass segrega- 
tion in these clusters again reduces tec, and for an 
IMF with a slope of —2.1, as may be p resent in 
the Arches cluster ( Espinoza et al.l 20091 ) , this re- 
duced tec is nearly the same as for an unsegregated 
cluster with a Salpeter IMF. Thus, if the IMF 
is Salpeter like, primordial mass-segregation in- 
creases the parameter space for runaway collisions 
to occur, whereas for flatter IMFs the parameter 
space remains very similar to that in unsegregated 
clusters with a Salpeter like IMF. From this we can 
conclude that if young massive star clusters (like 
Wester lund 1) are primordially mass segregated 
and have an IMF slope similar to Arches then the 
results from Paper I are directly applicable. 

We thank John Fregeau and Sourav Chatter- 
jee for the useful discussions. This work was sup- 
ported by NSF Grant AST-060.7498 and NASA 
Grant NNX08AG66G at Northwestern University. 
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Fig. 3. — Enclosed stellar mass in primordially 
mass-segregated (solid lines) and unsegregated 
(dashed lines) clusters. The average mass in the 
cluster is « 0.68 M©. The black and red solid 



lines a re for mass-segregated clusters (jSubr et al 



(|2008l ) recipe) with the degree of mass segrega- 
tion, S = 0.2 and S = 0.3, respectively. We 
show here mass enclosed in stars more massive 
than average mass of the cluster (Af, > 20 M© and 
M» > 50 Mq) as well as mass enclosed in stars of 
mass less than average mass (M* < 0.4 M©). All 
the simulations shown here contain 5 x 10^ stars 
with a Salpeter IMF (within Mmin = 0.2 M© and 
Mmax = 120 M©). Lines denoting enclosed stellar 
mass contained in stars more massive than the av- 
erage mass (i.e., > 20 M© and M, > 50 M©) 
show that the massive stars in the unsegregated 
clusters are at significantly larger radii than the 
primordially mass-segregated cluster. The me- 
dian position of massive stars in the unsegregated 
cluster differ by a factor of '-^ 4 for stars with 
A/» > 50 M©, and by a factor of ~ 3 for stars with 
M* > 20 M©, compared to the mass-segregated 
cluster with S — 0.3. Whereas lines denoting en- 
closed stellar mass contained in stars less massive 
than the average mass (i.e., < 0.4 M©), show 
that the low-mass stars are definitely much more 
spread out in case of primordially mass-segregated 
clusters when compared to the unsegregated clus- 
ter. 
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Fig. 4. — Same as FigjUb ut for primordially mass - 
segregated clusters using IBaumgardt et al.l (l2008l) 
recipe. The black solid lines denote the primor- 
dially mass-segregated clusters and the dashed 
lines are for unsegregated clusters. Massive stars 
in unsegregated clusters are at significantly larger 
radii, with their median position differing by a fac- 
tor of ^ 4 for stars with M^, > 50 M©, and by a fac- 
tor of ~ 3 for stars with > 20 M©, compared 
to the primordially mass-segregated clusters. 
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Fig. 5. — Parameter space survey for the pri- 
mordially unsegregated clusters. All the simu- 
lations shown here contain 5 x 10^ stars with a 
Salpeter IMF (within Mi„i„ = 0.2 Mq and M,„ax = 
120 Mq). Each circle represents a simulation of a 
cluster with a particular virial radius (rvir) and 
number of stars (N). Filled circles represent sim- 
ulations that result in a runaway while the empty 
circles indicate cases in which no runaway occur. 
The straight line denotes tec = 0.15irc = 3Myr 
as a function of rvir and N; clusters with initial 
conditions below this line have t^c < 3 Myr. As 
expected, we see in this plot that all the simula- 
tions that result in a runaway do in fact fall below 
this line. 



Fig. 6. — Same as FigE] but all the simula- 
tions are for primordially mass-segregated clusters 
{S = 0.3). The solid straight line once again de- 
notes denotes tec = 0.15irc = 3 Myr as a func- 
tion of Tvir and N for initially unsegregated clus- 
ters. This plot show that simulations of primor- 
dially mass-segregated clusters with initial condi- 
tions even above the solid straight line, might have 
core collapse time less than 3 Myr and can result 
in a runaway. Thus, simulations of primordially 
mass-segregated clusters increase the parameter 
space (range of rvir in pc) for runaway collisions to 
happen when compared to simulations of primor- 
dially unsegregated clusters by at least a factor of 
« 3. 
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Fig. 8. — Decrease in the core collapse time (tec) 
of a cluster with = 5 x 10^ stars, as the degree 
of primordial mass segregation (S) is increased in 
the cluster. For each value of S, the three different 
colors represent simulations with three different 
virial radii (rvh ). 



Fig. 7. — Growth curves of the most massive 
star (Mmms) of an unsegregated cluster (dashed 
lines) and primordially mass-segregated cluster 
(solid lines) with = 5 x 10^ stars and virial 
radii of 0.64 pc. The different colors indicate sim- 
ulations of the same cluster with different ran- 
dom seeds. While the most massive star in the 
primordially mass-segregated cluster have a very 
steep and rapid mass growth leading to a runaway, 
the most massive star in the unsegregated clus- 
ter shows no such effect. For the unsegregated 
cluster = 38.67(red), 38.33(blue), 38.66(green) 
and 40 Myr (purple) while for the primordially 
mass-segregated cluster tic values are 13.53(red), 
14.0(blue), 14.0(green) and 13.33 Myr(purple). 




Fig. 9. — Dependence of core collapse time in units 
of the central relaxation time, on primordial mass 
segregation (S*). For each value of S, the three 
different colors represent simulations with three 
different virial radii (rvir). For comparison, tec /tic 
as predicted in Paper I for unsegregated clusters is 
plotted as the dotted line. Even for the mass seg- 
regated clusters tcc/^ic remains roughly consistent 
with Paper I. 
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Fig. 10. — Decrease in the central relaxation time 
{tic in Myr) of a cluster of 5 x 10^ stars with the 
increase in the degree of primordial mass segrega- 
tion (S). For each value of S, the three different 
colors represent simulations with three different 
virial radii (rvh ). 
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Fig. 11. — Evolution of the central relaxation time 
tic (solid line) and the factors affecting t^c (dot- 
ted lines). All parameters are normalized to their 
values for an unsegregated cluster, (X/Xg=o.o)- 
The black line is for the inverse of central den- 
sity (1/ pc), the intermediate grey is for the inverse 
of average stellar mass in the core (1/ {m)^), and 
the light grey represents the cube of the velocity 
dispersion in the core (cr^)- While remains al- 
most the same with increasing S, both 1/pc and 
1/ {m) ^ decrease ^rc- However, for higher values of 
S' (5' > 0.15), it is mainly 1/ (m)^ which decreases 
iic • 



0.003 




0.05 0.1 0.15 0.2 0.25 0.3 



Fig. 12. — Fraction of the total stellar mass in the 
cluster, that ends up as a VMS after a runaway 
collision, as a function of the degree of mass seg- 
regation. All the simulations are for a cluster of 
5 X 10^ stars and a virial radius of 0.64 pc. The 
■^VMS /Mtot shown here is averaged over 5 simula- 
tions with the same initial conditions and different 
random seeds. The error bars are the standard de- 
viations calculated from the data and the average 
value. 
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Fig. 13. — Dependence of the core collapse time 
in units of the half-mass relaxation time, on the 
shape and width of the IMF. The horizontal axis 
shows the ratio of maximum to mean stellar mass 
{mmaxi (w)) in the IMF. Th e solid black line 



i s for mass-segregated clusters (jBaumgardt et al 



whereas all other symbols are for unsegre- 
gated clusters. The black dashed straight line, the 
triangles and the stars denote the simulations with 
physical stellar collisions turned off whereas the 
circles represent simulations with collisions (trun- 
cated Salpeter IMF with M^ax = HO, 100, 90, 80, 
70, 60). The stars correspond to the Variable IMF 
varied as a power law from —1.72 to —2.3 {aA) and 
the triangles correspond to the middle section of 
the Variable IMF varied as a power law from —1.9 
to —2.3 (a3). Flatter IMFs (triangles and stars) 
increase ice by 30% but when collisions (circles) 
are included t^c decrease on an average by 15%. 
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Table 1 
CLASSIFICATION OF IMF 



IMF 


Range of m 


a 








Kroupa 








0.1-0.5 


-1.3 




0.5 - oo 


-2.3 


Arches 








0.1-0.5 


-1.3 




0.5 - 1.0 


-2.3 




1.0-3.0 


-2.04 




3.0-15.0 


-1.5 




15.0 -00 


-1.72 


Variable 








0.1-0.5 


-1.3 




0.5 - 1.0 


-2.3 




1.0-3.0 


-2.04 




3.0 - CXD 


-1.72 



21 



Table 2 

INITIAL CONDITIONS AND RESULTS FOR SIMULATIONS WITH VARYING IMF. 



Name 




IMF 




^vir 

{pc) 




01 


Variable 


0:3 = - 


-2.04,04 = 


-1.7 


0.39 


0.127 


02 


Variable 


as = 


-2.04,0:4 = 


-1.8 


0.39 


0.103 


03 


Variable 


as = 


-2.04,04 = 


-1.9 


0.39 


0.092 


04 


Variable 


as ^ - 


-2.04,04 = 


-2.0 


0.39 


0.088 


05 


Variable 


as = - 


-2.04,04 = 


-2.1 


0.39 


0.077 


06 


Variable 


as = - 


-2.04,04 = 


-2.2 


0.39 


0.081 


07 


Variable 


a^ = 


-2.04,04 = 


-2.3 


0.39 


0.078 


08 


Variable 


as = 


-1.9, 04 = - 


-1.72 


0.39 


0.132 


09 


Variable 


as = - 


-2.0,04 = - 


-1.72 


0.39 


0.125 


10 


Variable 


as = - 


-2.1, 04 = - 


-1.72 


0.39 


0.117 


11 


Variable 


as = 


—2.2, 04 = - 


-1.72 


0.39 


0.102 


12 


Variable 


as = 


-2.3, 04 = - 


-1.72 


0.39 


0.099 


13 


Salpeter 




= 110 




0.39 


0.067 


14 


Salpeter 




= 100 




0.39 


0.073 


15 


Salpeter 




= 90 




0.39 


0.0812 


16 


Salpeter 


M„,ax 


= 80 




0.39 


0.079 


17 


Salpeter 




= 70 




0.39 


0.097 


18 


Salpeter 


Mmax 


= 60 




0.39 


0.108 



'^All models start as Plummer spheres with = 5 x 10^ 
and the Variable IMF covers the same mass range (0.2 — 
120Mq). 



Note. — rvir is the virial radius of the cluster, tec denotes 
the core collapse time of the cluster and trh denotes the half 
mass relaxation time. 
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Table 3 

INITIAL CONDITIONS AND RESULTS OF SIMULATIONS WITH PRIMORDIAL MASS 

SEGREGATION (5UBR RECIPE). 



Name 


N 


^vir 


S 


irc(O) 


^CC 






(pc) 




(Myr) 


(Myr) 


la 


3x10^ 


0.50 


0.00 


29.30 


4.40 


lb 


3x10^ 


0.50 


0.05 


21.86 


3.28 


Ic 


3x10^ 


0.50 


0.10 


19.33 


2.90 


Id 


3x10^ 


0.50 


0.20 


11.26 


1.69 


le 


3x10^ 


0.50 


0.30 


6.80 


1.02 


2a 


3 xlO^ 


0.62 


0.00 


28.66 


4.30 


2b 


3x10^ 


0.62 


0.05 


24.53 


3.68 


2c 


3x10^ 


0.62 


0.10 


13.00 


1.95 


2d 


3x10^ 


0.62 


0.20 


10.33 


1.55 


2e 


3x10^ 


0.62 


0.30 


5.86 


0.88 


3a 


3 xlO^ 


0.77 


0.00 


41.33 


6.20 


3b 


3x10^ 


0.77 


0.05 


31.53 


4.73 


3c 


3x10^ 


0.77 


0.10 


20.80 


3.12 


3d 


3x10^ 


0.77 


0.20 


19.20 


2.88 


3fe 


3x10^ 


0.77 


0.30 


13.93 


2.09 


5a 


5x10^ 


0.41 


0.00 


25.86 


3.88 


5b 


5x10^ 


0.41 


0.05 


21.13 


3.17 


5c 


5x10^ 


0.41 


0.10 


17.4 


2.61 


5d 


5x10^ 


0.41 


0.15 


10.00 


1.50 


5e 


5x10^ 


0.41 


0.20 


8.66 


1.30 


5f 


5x10^ 


0.41 


0.25 


8.0 


1.20 


5g 


5x10^ 


0.41 


0.30 


5.53 


0.83 


6a 


5x10^ 


0.52 


0.00 


31.86 


4.78 


6b 


5x10^ 


0.52 


0.05 


22.53 


3.38 


6c 


5x10^ 


0.52 


0.10 


13.93 


2.09 


6d 


5x10^ 


0.52 


0.15 


14.13 


2.12 


6e 


5x10^ 


0.52 


0.20 


10.46 


1.57 


6f 


5x10^ 


0.52 


0.25 


10.53 


1.58 


6g 


SxlO'"^ 


0.52 


0.30 


9.93 


1.49 


7a 


5x10^ 


0.64 


0.00 


38.66 


5.80 


7b 


5x10^ 


0.64 


0.05 


23.33 


3.50 


7c 


SxlO'"^ 


0.64 


0.10 


16.66 


2.50 


7d 


5x10^ 


0.64 


0.15 


17.53 


2.63 


7e 


5x10^ 


0.64 


0.20 


16.00 


2.40 


7f 


SxlO'"^ 


0.64 


0.25 


16.00 


2.40 


7g 


5x10^ 


0.64 


0.30 


13.53 


2.03 


7h 


5x10^ 


0.64 


0.00 


38.33 


5.75 


7i 


5x10^ 


0.64 


0.05 


23.33 


3.50 


7j 


5x10^ 


0.64 


0.10 


15.93 


2.39 


7k 


5x10^ 


0.64 


0.15 


17.86 


2.68 
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Table 3 — Continued 



AT 

Name 


TV 




s 




^CC 






(pc) 




(Myr) 


(Myr) 


71 


5x10^ 


0.64 


0.20 


16.00 


2.40 


7m 


5x10^ 


0.64 


0.25 


15.93 


2.39 


7n 


5x10^ 


0.64 


0.30 


14.00 


2.10 


7o 


5x10^ 


0.64 


0.00 


38.66 


5.80 


7d 


5x10^ 


0.64 


0.05 


25.33 


3.80 




5x10-'' 


0.64 


0.10 


16.66 


2.50 


7r 


5x10^ 


0.64 


0.15 


10.76 


2.69 


7s 


5x10^ 


0.64 


0.20 


16.00 


2.40 


7t 


5x10-'' 


0.64 


0.25 


16.00 


2.40 


7u 


5x10^ 


0.64 


0.30 


14.06 


2.11 


7v 


5x10^ 


0.64 


0.00 


40.00 


6.00 


7x 


5x10''' 


0.64 


0.05 


23.33 


3.50 


T'w 


5x10^ 


0.64 


0.10 


16.86 


2.53 


7v 


5x10^ 


0.64 


0.15 


17.53 


2.63 


7zl 


5x10''' 


0.64 


0.20 


16.00 


2.40 


7z2 


5x10^ 


0.64 


0.25 


16.00 


2.40 


7z3 


5x10^ 


0.64 


0.30 


13.33 


2.00 


8a 


5x10^ 


0.77 


0.00 


52.60 


7.89 


8b 


5x10^ 


0.77 


0.05 


34.60 


5.19 


8c 


5x10^ 


0.77 


0.10 


28.26 


4.24 


8d 


5x10^ 


0.77 


0.15 


24.20 


3.63 


8e 


5x10^ 


0.77 


0.20 


20.40 


3.06 


8f 


5x10^ 


0.77 


0.25 


19.73 


2.96 


8e- 


5x10''' 


0.77 


0.30 


16.53 


2.48 


9a 


1x10^ 


0.41 


0.10 


20.8 


3.12 


9b 


1x10^ 


0.41 


0.15 


19.06 


2.86 




1 y1 

J_XXU 


U .4: 1 


u.zu 


±u.oo 




9d 


IxlQS 


0.41 


0.30 


10.40 


1.56 


10a 


IxIQS 


0.51 


0.10 


32.60 


4.89 


10b 


IxlQf^ 


0.51 


0.15 


21.46 


3.22 


10c 


IxlQS 


0.51 


0.20 


17.66 


2.65 


lOd 


1x10*5 


0.51 


0.30 


15.53 


2.33 


11a 


IxlQS 


0.83 


0.10 


46.53 


6.98 


lib 


IxlQS 


0.83 


0.15 


43.73 


6.56 


He 


1x10*5 


0.83 


0.20 


27.60 


4.24 


lid 


1x10*^ 


0.83 


0.30 


20.06 


3.01 


12a 


1x10** 


1.10 


0.10 


60.53 


9.08 


12b 


1x10*5 


1.10 


0.15 


59.26 


8.89 


12c 


1x10** 


1.10 


0.20 


52.40 


7.86 


12d 


1x10** 


1.10 


0.30 


44.46 


6.67 
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Table 3 — Continued 



Name Tvir S irc(O) ice 

(pc) (Myr) (Myr) 



Note. — Here N and rvir are the initial 
number of stars and the initial virial radius 
respectively. S is the degree of mass seg- 
regation (equation ([SJ), tic(O) is the initial 
central relaxation time and tec denotes the 
core collapse time. 

^All initial models are isolated Plummer 
spheres with standard Salpeter IMFs with 
mass range (0.2 - I2OM0). 

'^Calculation of tjc and t^c- For Plum- 
mer models — 1 and = 0.769 
in A^-body units (Paper I, their Table 
1). From Tv in physical units we calcu- 
late in physical units. With initial N 
and initial (in physical units) we also 
calculate the Initial half-mass relaxation 
time in physical units, given by t^h = 
{0.138N/ln-fN)x{rl/GM)^/^. From U 
we calculate t^c {trh and ire for Plummer 
models is 0.093 and 0.0437 respectively, in 
Fokker-Planck units). Finally we obtain 

^cc ^ 0.15tic- 
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